Density-functional calculation of ionization energies of current-carrying atomic states 



m 
O 
O 

(N 

G 
G 



E. Orestes, T. Marcasso, and K. Capelle* 
Departamento de Qmmica e Fisica Molecular 
Institute: de Qmmica de Sao Carlos 
Universidade de Sao Paulo 
Caixa Postal 780, Sao Carlos, 13560-970 SP, Brazil 
(Dated: February 2, 2008) 

Current-density-functional theory is used to calculate ionization energies of current-carrying 
atomic states. A perturbative approximation to full current-density-functional theory is imple- 
mented for the first time, and found to be numerically feasible. Different parametrizations for the 
current-dependence of the density functional are critically compared. Orbital currents in open-shell 
atoms turn out to produce a small shift in the ionization energies. We find that modern density 
functionals have reached an accuracy at which small current-related terms appearing in open-shell 
configurations are not negligible anymore compared to the remaining difference to experiment. 
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I. INTRODUCTION 

Density-functional theory (DFT) [1] is based on identi- 
fication of the single-body charge density n(r) as key vari- 
able in terms of which all ground-state observables of an 
interacting many-electron system can be described [2, 3]. 
Although n(r) is in principle sufficient for this purpose, 
in practice it often turns out to be advantageous to em- 
ploy additional variables. The most commonly employed 
such additional variable is the spin magnetization m(r) 
(or the spin-resolved charge densities n-f(r) and n^(r)), 
leading to spin-density functional theory (SDFT) [2, 3, 4]. 
Other choices have occasionally been found useful, e.g. in 
solid-state physics [5, 6]. 

In the present paper we are interested in orbital mag- 
netism produced by currents forming in open-shell atoms. 
The current density seems a useful extra variable in this 
situation. In the absence of external magnetic fields 
ground-state currents can in principle be calculated also 
by means of (S)DFT. However, in practice the calcula- 
tion of orbital currents and their effects on observables is 
impossible in standard formulations of DFT and SDFT, 
because no explicit prescription for calculating the true 
(many-body) current density is known in these theories. 
Of course, one can always calculate the current aris- 
ing from the single-particle Kohn-Sham (KS) orbitals of 
SDFT or DFT, but these orbitals are constrained only to 
reproduce the correct charge (and spin) densities of the 
interacting many-body system, and there is no guarantee 
that the current calculated from them bears any relation 
to the true current. 

Orbital magnetism is thus basically out of reach of con- 
ventional DFT. In view of the importance of currents, 
appearing cither spontaneously or induced by external 
magnetic fields in a wide variety of many-body systems, 
such as atoms and molecules with open shells, nuclear- 



magnetic resonance, cyclotron resonance, superconduc- 
tivity, and magnetism of transition metal and rare-earth 
compounds, it is thus desirable to develop a DFT-based 
approach that allows to directly address the effects of 
orbital currents. In the present paper we explore one 
such formulation of DFT, namely nonrelativistic current- 
density functional theory (CDFT) [7, 8, 9]. Relativistic 
DFT [10, 11, 12] in principle also provides explicit infor- 
mation on the current, but standard implementations of 
it are formulated in a spin-only version, which prohibits 
extraction of information on the currents. Furthermore, 
the formalism of relativistic DFT is considerably more 
complicated than that of SDFT. 

CDFT is formulated in terms of the charge density 
n(r) and the nonrelativistic paramagnetic current density 
vector 

j P (r) = ^2i P ,k(r) 

k 

= 2^? E MktoV^r) " (V^.(r))^(r)] . (1) 



This current is to be calculated from the CDFT Kohn- 
Sham equations [7, 8, 9] 



where an upper index 'c' denotes CDFT, 

,2 



V s %r) = v c s (r) 



2mc" 



{A ext {vf - A s {vf 



and 



<(r) = v ext {r) + v H (r) + v c xc (r) 



A s (r) = A ext (r) + A xc (r). 



(2) 



(3) 



(4) 



(5) 
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Here v ex t and A ex ± are external static electric and mag- 
netic potentials, vh is, as usual, the Hartree potential, 
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and v xc and A xc are the exchange-correlation (xc) scalar 
and vector potentals of CDFT, respectively [7, 8, 9]. 
Gauge invariance of CDFT requires that the xc energy 
E X c[n,3p} depends on the current only through the so 
called vorticity [7, 8, 9] 

*jp]W=Vx J ^, (6) 

i.e., is of the form E xc [n,j p ] — E xc [n,v\. This depen- 
dence provides a useful constraint on approximate CDFT 
functionals. 

Although CDFT formally solves the problem how 
to obtain current-related information from DFT, many 
practical questions remain. One is, clearly, the construc- 
tion of approximate current-density functionals. A brief 
summary of progress in this area is given in Sec. II. An- 
other is the actual implementation of CDFT. In practice, 
a fully self-consistent solution of the CDFT equations is 
still quite demanding. Furthermore, in many interest- 
ing situations the effect of orbital magnetism, while im- 
portant for a qualitative and quantitative understanding, 
is relatively small, so that a fully self-consistent treat- 
ment of the orbital degrees of freedom is not always re- 
quired. The question thus arises if one can put the in- 
sights, gained from the formalism of full-fledged CDFT, 
at work within conventional DFT, in order to achieve 
an improved description of orbital magnetism via a self- 
consistent solution of the widely implemented traditional 
Kohn-Sham equation. 

A simple answer to this question was given by one of us 
in Ref. [13], by pointing out that the CDFT Kohn-Sham 
equations can be written in the form of the DFT ones 
plus a remainder that depends explicitly on the xc and 
external vector potentials. In the absence of external 
magnetic fields and for not too large xc vector poten- 
tials it then suggests itself to use low-order perturbation 
theory in order to describe the CDFT modifications to 
the DFT equations. This idea has been worked out in 
Ref. [13], where explicit expressions for the CDFT cor- 
rections to a number of important DFT quantities were 
given. The resulting approach is labeled 'perturbative 
CDFT', or simply pCDFT. A simple example of pCDFT 
expressions is the correction to the DFT eigenvalue spec- 
trum ef, which can be cast in the form [13] 

4 =4~ \ J d 3 r££(r) ■ A xc [n,$ s ](r). (7) 

II. CURRENT-DEPENDENT LOCAL-DENSITY 
APPROXIMATION 

In the present paper we apply pCDFT to a study of 
the effect of currents in open shells on atomic ionization 
energies. To this end, we require an explicit expression 
for the xc energy in the presence of orbital currents. In 
Rcfs. [7, 8] Vignale and Rasolt proposed an extension 
of the local-density approximation (LDA) of ordinary, 



charge-only, DFT to the case of CDFT. Their functional 
takes the form 

E xc [n,j p ] = J d 3 rn(r)e xc (n,u) , (8) 

where 

exc(n,i/) = e IC (n,0) + ^^-ljiA^-. (9) 

Here k F is the Fermi wave vector of a non-interacting 
electron gas, related to the density via k F = (Z^n) 1 ^ . 
The function e xc (n, 0) is the exchange-correlation energy 
per particle in the absence of external magnetic fields, 
and can be approximated, e.g., by the usual LDA or 
any of the available generalized-gradient approximations 
(GGA's). The second term on the right-hand side of 
Eq. (9) is a weak-field (linear response) expression for 
the current-dependent part of the functional. The func- 
tional defined by the preceeding equations has been used 
in several CDFT calculations [14, 15, 16, 17] and should 
be good enough for a first orientation about the size and 
nature of current-related phenomena. 

Many-body effects enter this functional via the ratio of 
the orbital susceptibilities of the interacting and the non- 
interacting electron gas, s :— Xl/Xl- This ratio has been 
calculated numerically by Vignale, Rasolt and Geldart 
[18] for 10 values of the dimensionless density parameter 
r s , which is related to the density by n — 3/47r(r s ao) 3 , 
where ao is the Bohr radius. To utilize these results in 
CDFT, Lee, Colwell and Handy (LCH) proposed the ex- 
pression [16, 17] 

s L CH{r s ) = (1.0 + 0.028r s ) exp (-0.042r s ) (10) 

as a convenient and accurate analytical interpolation 
through the numerical data of Ref. [18]. In the origi- 
nal reference [16] this fit is claimed to have an rms error 
of ~ 1.5 x 10~ 3 . This appears to be a misprint in [16], be- 
cause on redoing the calculation we find that over the 10 
data points provided by Vignale, Rasolt and Geldart [18] 
the rms error of expression (10) is a little bigger, namely 
~ 1.9 x 10~ 3 . We have also constructed two alternative 
interpolations through the data of Ref. [18]. With the 
same number of fitting parameters (three) as in the LCH 
fit we find that the expression 

s 3 {r s ) = 0.9956 - 0.01254r s - 0.0002955^ (11) 

has an rms error of ~ 1.2 x 10~ 3 , while the 5-term ex- 
pression 

s 5 {r s ) = 1.1038 - 0.4990r s 1/3 + 0.4423^ 

-0.06696r s + 0.0008432^ (12) 

has an rms error of only <~ 2.1 x 10~ 4 over the same 10 
data points. Over the range of values of r s spanned by 
these data (1 . . . 10) the latter expression should thus be 
prefered, compared to (10) or (11). In the calculations 
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presented below the current flows mainly in a region in 
which 0.1 < r s < 5. We have therefore performed all 
calculations once with the LCH fit (10) and once with 
the above 5-term fit (12). As will be seen below, the 
differences between both sets of results can be consider- 
able. Since our calculations are self-consistent only with 
respect to the charge density, but perturbative with re- 
spect to the current density, we expect, however, that it 
is only the order of magnitude of the results that is quan- 
titatively reliable, and on this order of magnitude both 
parametrizations employed agree consistently. 

The LCH expression and both alternative fits differ 
markedly for r s > 12, i.e., in the extreme low-density 
limit. For the present calculation this range is less im- 
portant, but in view of potential future applications of 
CDFT in the low-density regime it should be noted that 
the numerical data of Ref. [18] do not constrain the vari- 
ous fits in that regime. In the opposite limit, r s — > 0, the 
five-term fit (12), although on average significantly more 
accurate than the LCH expression, does not correctly re- 
cover the value at r s = 0, whereas the LCH fit correctly 
yields slch(0) = 1- However, the limit r s — > corre- 
sponds to infinite density, and is thus rather unimportant 
for typical atomic physics applications. Moreover, in this 
limit the asymptotically exact expression [18] 

s re ^ (r s ) = 1 + —r s lnr s + 0.08483— + 0(r 2 a ), (13) 

D7T 7T 

where a = (4/97T) 1 / 3 , is available, so that there is no need 
for data fitting at all in this limit. 

In addition to Eq. (8) with (9) several other CDFT 
functionals have been proposed, but none seems suitable 
for our purposes. The approach of Ref. [21] has not yet 
led to an explicit expression for the current dependence 
of A xc [n, jp]. The functional of Ref. [22] displays quan- 
tum oscillations arising from Landau-level filling in the 
electron gas, and is thus suitable only for extended sys- 
tems. The functionals of Refs. [23, 24], on the other hand, 
were designed specifically for two-dimensional systems in 
the quantum Hall regime. The expression of Ref. [25], 
finally, is not a vorticity functional. We will return to 
this last functional in our discussion of the results in Ta- 
ble II, below. Given the scarcity of suitable functionals, 
the proposal of Ref. [26] to generate a CDFT functional 
by means of a set of integral transformations from an in- 
put SDFT functional may prove useful in the future. In 
the present work, however, we restrict attention to the 
simple linear- response LDA defined by Eqs. (8) and (9). 



III. NUMERICAL RESULTS FOR IONIZATION 
ENERGIES OF CURRENT-CARRYING ATOMIC 
STATES 

After this preparatory discussion of the CDFT func- 
tional, we now return to the question of CDFT shifts 
with respect to the DFT eigenvalues. In general the 
eigenvalues obtained from solution of the KS equations 



have no rigorous physical meaning, although they can 
bear a semiquantitative relationship to the true energy 
spectrum [27]. One exception to this rule is the highest 
occupied eigenvalue of the KS spectrum, which is known 
to be the negative of the system's ionization energy, and 
as such can be compared directly to experiment. Al- 
though once the system is prepared in a current-carrying 
state there are, in general, CDFT shifts to all KS eigen- 
values, we therefore focus in the present paper on the 
highest occupied one, since it is most tightly connected 
with experiment. 

In comparisons of experimental ionization energies 
with KS eigenvalues obtained within the local-density ap- 
proximation (LDA) it typically turns out that the LDA 
eigenvalues are significantly off. The origin of this prob- 
lem is known to be the asymptotic behaviour of the LDA 
xc potential, which decays too rapidly, and leads to too 
weak binding of the outermost electron. Self-interaction 
corrections (SIC) [28], which correct the wrong LDA 
asymptotics, significantly improve on the LDA values for 
ionization energies, and typically are quite close to ex- 
perimental data [19]. Once a self-interaction correction 
has been applied, and the main error of the LDA ion- 
ization energies removed, one can consider the effect of 
additional small corrections, such as the effects of orbital 
currents. Such current-dependent corrections to the ion- 
ization energy can become important, e.g., when ioniza- 
tion (or transfer of electrons during formation of chemical 
bonds) takes place in the presence of external magnetic 
fields, since such fields polarize the atom and can give 
rise to orbital currents. 

It is in this situation, ionization or electron transfer 
in the presence of static external magnetic fields, where 
pCDFT calculations of current-induced shifts of the ion- 
ization energies are directly applicable. Within the con- 
text of density-functional theory, the possibility of small 
current-induced shifts in the ionization energies is also 
relevant for the calculation of excited states from time- 
dependent DFT (TDDFT). A recent systematic investi- 
gation [29] of sources of error in excitation energies cal- 
culated from TDDFT concludes that the most impor- 
tant ingredient in such calculations is the ground-state 
xc potential used to generate the Kohn-Sham response 
function, and that 'the most important requirement for 
such a potential would be that its highest occupied eigen- 
value reproduces the experimental ionization potential as 
closely as possible' (cf. Sec. 6 of Ref. [29]). TDDFT 
calculations of excitations from open shells are thus ex- 
pected to sensitively depend on current-induced shifts of 
the ionization energies, if the excitation takes place in 
the presence of external magnetic fields. 

A separate issue is whether a CDFT calculation of the 
type presented here can also be useful in the absence of 
any external magnetic field. An argument for such utility 
could run along the following lines: The negative of the 
highest occupied KS eigenvalue of DFT gives the true 
ionization energy. The negative of the highest occupied 
eigenvalue of CDFT also gives this energy. In the ab- 
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TABLE I: Current-induced changes in the ionization energies of atoms with open p (B to CI) and d (Sc, Y) shells. First column: 
atom. Second column: zero-current ionization energy calculated within LDA-SIC-KLI, from [19]. Third column: experimental 
ionization energies, from [20]. Fourth column: current-carrying single-particle state considered here. Fifth column: pCDFT 
correction in LDA, using the LCH expression (10) for the susceptibilities. Sixth column: pCDFT correction in LDA, using the 
present expression (12) for the susceptibilities. All values in eV. 





I 

LDA-SIC-KLI 


I 

exptl. 


cc sp state(s) 




-Ae pCDFT 
with (10) 


_^ e pCDFT 

with (12) 


B 


8.316 


8.2980 


m = 1 




0.072 


0.056 


C 


11.60 


11.2603 


mi = 1, ni2 = 




0.045 


0.051 








mi = 1, 7712 = 1 




0.18 


0.20 


N 


14.95 


14.5341 


mi = 1, 7712 = 0, 7773 


= 


0.034 


0.047 








mi = 1, 7772 = 1, 7773 


= 


0.14 


0.19 





14.33 


13.6181 


mi = 1, 7772 = 1, 7773 


= 0,m 4 = 


0.11 


0.11 








mi = 1, 7712 = 1, 7773 


= 0, 7774 = —1 


0.027 


0.027 


F 


18.61 


17.4228 


mi = 1, 7772 = 1, 7773 


= 0, 7774 = 0, 7774 = —1 


0.023 


0.040 


Al 


5.570 


5.9858 


m = 1 




0.049 


0.030 


Si 


7.804 


8.1517 


mi = 1, 7712 = 




0.022 


0.018 








mi = 1, 7712 = 1 




0.089 


0.073 


P 


10.07 


10.4867 


mi = 1, 7712 = 0, 7773 


= 


0.015 


0.014 








mi = 1, 7712 = 1, 7773 


= 


0.059 


0.057 


S 


10.41 


10.3600 


mi = 1, 7712 = 1, 7773 


= 0, 7774 = 


0.044 


0.044 








mi = 1, 7712 = 1, 7773 


= 0, 7774 = — 1 


0.011 


0.011 


CI 


13.08 


12.9676 


mi = 1, 7772 = 1, 7773 


= 0, 7774 = 0, 7774 = —1 


0.0091 


0.011 


Sc 




6.5615 


m = 1 




0.036 


0.043 








m = 2 




0.037 


0.044 


Y 




6.2171 


m = 1 




0.035 


0.035 








m = 2 




0.036 


0.037 



sence of external magnetic fields both values must thus 
be identical, if one works with the exact functionals. In 
practise, of course, we do not have the exact functionals 
available, and it becomes a meaningful question to ask 
which of the two approximate eigenvalues is closer to the 
experimental energy. Our original expectation (not con- 
sistently confirmed by the numerical data shown below) 
was that for atoms whose many-body ground state has 
a nonzero value of the total angular momentum quan- 
tum number L the CDFT value would be better than 
the DFT one, if the CDFT-KS calculation is performed 
for a system prepared in a state with an orbital current. 

Motivated by these considerations we have numerically 
calculated the pCDFT correction to the highest occupied 
KS eigenvalue for a series of atoms with unfilled p and d 
shells, prepared in a current-carrying state. Energies and 
radial wave functions of the unperturbed system were 
obtained numerically from a standard spherically aver- 
aged DFT calculation, using the basis-set-free program 
opmks [31], both within the LDA and the BLYP GGA 
[32]. The full single-particle orbital for a current-carrying 
state was then obtained by multiplying the radial wave 
function with a spherical harmonic corresponding to a 
definite value of the magnetic quantum number m. Such 
states carry a current proportional to m. This proce- 



dure defines how we prepare a current-carrying state in 
the KS system. Experimentally, the preparation of an 
open-shell atom in a current-carrying state is achieved 
by applying suitable external magnetic or electric fields, 
as in the Zeeman and Stark effects [33] . 

In Table I the results of LDA-SIC and pCDFT calcula- 
tions of atomic ionization energies are listed for a number 
of current-carrying states. It is interesting to see that the 
differences between values for Ae pCDFT obtained with 
both parametrizations of the susceptibility data used in 
this work, Eq. (10) and Eq. (12), can be substantial, al- 
though the order of magnitude predicted by both is the 
same. On the other hand, it makes little difference for the 
size of Ae pCDFT whether the unperturbed orbitals used 
for calculating j Pi/ t(r) and A xc (r) are obtained from LDA 
or GGA. The values listed in Table I were obtained with 
LDA. For comparison we also performed calculations us- 
ing the BLYP GGA [32] , but the resulting changes in the 
pCDFT corrections are consistently smaller than the ones 
arising from changing the paramctrization of the current- 
dependent part of the functional from (10) to (12). On 
the other hand, depending on the particular occupation 
of m-substates in the noninteracting system the size of 
both the current and the resulting energy shift can vary 
considerably. To explore this variation Table I lists values 
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for more than one occupation for most atoms. 

Comparison of the experimental ionization energies 
listed in Table I with the sum of LDA-SIC-KLI values 
and the pCDFT correction shows that pCDFT slightly 
improves agreement with experiment for several second- 
row elements, but worsens it for the first-row ones. The 
data are thus inconclusive as to whether pCDFT repre- 
sents a true improvement on calculations of ionization 
energies that neglect currents. This is similar to the sit- 
uation encountered by LCH in the CDFT calculation of 
nuclear shielding tensors, in Ref. [17], where the same 
functional (9) as here was employed and CDFT was found 
to yield rather small shifts that did not systematically 
improve agreement with experiment. 

In Ref. [17] this was attributed to deficiencies in the 
approximation for the density functional. Although the 
same could be true here, too, it is in our opinion more 
likely that the inconclusive comparison between LDA- 
SIC-KLI+pCDFT ionization energies with experiment is 
due to the fact that while the pCDFT values were ob- 
tained by explicitly assuming a current-carrying config- 
uration in the open shell, the experimental data refer to 
configurations with on average equal population of all 
(2L + 1)(2S + 1) degenerate substates belonging to the 
2S+1 L ground-state term of the atom. Such an average 
eliminates the orbital currents. A direct comparison of 
the sum of the LDA-SIC-KLI values and the pCDFT ones 
with available experimental data is thus not necessarily 
meaningful. Rather, experiments in the presence of suit- 
able external fields, selectively populating states with a 
nonzero orbital current, are called for. 

Even in the absence of such experiments, it is, how- 
ever, still possible to compare the size of the pCDFT 
shifts in Table I with other calculations of the effects of 
orbital currents in open-shell atoms and with the accu- 
racy of present-day DFT calculations. Such a comparison 
is presented in Table II, in which we compare the order 
of magnitude of the pCDFT shifts in Table I with the 
deviation of LDA-SIC calculations from experiment, and 
with results from another recent calculation of effects of 
orbital currents in open-shell atoms. From the numbers 
in Table II it is obvious that once the effect of the self- 
interaction corrections has been taken into account, the 
order of magnitude of the effect of currents in the open 
shell is comparable to the remaining difference between 
the LDA-SIC data and experiment. As explained above, 
this does not mean that agreement with experiment is 
necessarily improved by simply adding the two contribu- 
tions, but it implies that further refinement of density 
functionals for the calculation of the electronic structure 
of atoms must take the possibility of curent-dependcnt 
energy shifts into account [25, 29, 30]. 

In recent work of Becke [25] the effect of orbital cur- 
rents in open shells is studied from a different point 
of view. Instead of explicitly preparing the atom in 
a symmetry-broken current-carrying state and studying 
the resulting shifts of single-body energies, as we do here, 
he considers the total ground-state energies calculated 



TABLE II: Comparison of order of magnitude of current- 
related effects in open shells with accuracy of LDA-SIC cal- 
culations. First column: atom considered. Second column: 
absolute difference between LDA-SIC-KLI data [19] and ex- 
perimental result [20]. Third column: current- induced shift 
of ground-state energies, as calculated from fourth and fifth 
columns of Table 2 of Ref. [25]. Fourth column: current- 
related contribution to ionization energy of current-carrying 
states, taken from the last column of Table I. (For atoms for 
which more than one configuration is listed in Table I we have 
taken the one with the more negative Ae p for the present 

comparison.) Note that the numbers in columns two, three 
and four measure different things, and should not expected to 
be identical. As discussed in the main text, the fact to note 
is that they are of similar magnitude, showing that current- 
related phenomena in open-shell atoms (columns three and 
four) lead to energy shifts that are comparable to the preci- 
sion of modern density functionals (column two). All values 
in eV. 





deviation of 


-A£(j) 


_ Ae pCDFT 




LDA-SIC-KLI [19] 


from Ref. [25] 


with (12) 




from experiment [20] 






B 


0.018 


0.14 


0.056 


C 


0.34 


0.15 


0.20 


N 


0.42 




0.19 





0.71 


0.23 


0.11 


F 


1.2 


0.24 


0.040 


Al 


0.42 


0.075 


0.030 


Si 


0.35 


0.069 


0.073 


P 


0.42 




0.057 


S 


0.050 


0.092 


0.044 


CI 


0.11 


0.085 


0.011 



from different but symmetry-equivalent configurations, 
which should be degenerate. Nevertheless, approximate 
density functionals can give rise to an artificial break- 
ing of the degeneracy with respect to M [25, 34]. Once 
this degeneracy is broken, the corresponding ground state 
carries a current, just as in the present calculations. Ac- 
cording to [25] the original degeneracy is restored if the 
current dependence of the xc functional is explicitly ac- 
counted for. The current-dependent functional used in 
Ref. [25] is not constructed within the CDFT of Vig- 
nale and Rasolt but based on analysis of the dependence 
of the curvature of the xc hole on the Kohn-Sham cur- 
rent and, as pointed out above, is not a vorticity func- 
tional. Implementation of the functional is done not 
self-consistently (as in full CDFT) or perturbatively (as 
in pCDFT), but in a post-LDA manner, in which or- 
bitals obtained in a self-consistent (current-independent) 
LDA calculation are substituted once into the current- 
dependent functional [35] . In spite of these methodolog- 
ical differences between the present work and Ref. [25], 
Becke also finds that to within the accuracy of today's 
density-functionals the current-dependent terms cannot 
be neglected in high-precision DFT calculations of open- 
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shell atoms. In fact, as illustrated in the third column 
of Table II, the size of the current-dependent corrections 
to the total energy differences obtained from subtracting 
the values in the last two columns of Table 2 of Ref. [25] 
is rather similar to the size of the current-induced cor- 
rections to efc, calculated here. 

IV. CONCLUSIONS 

We have performed CDFT calculations of ionization 
energies of open-shell atoms prepared in a current- 
carrying state, with the aim to illustrate the usefulness 
and viability of CDFT for electronic-structure calcula- 
tions in the presence of orbital currents. We summarize 
our conclusions as follows: 

(i) The perturbative approximation to CDFT, pCDFT, 
has been implemented for the first time and was found 
to be numerically feasible. 

(ii) The CDFT xc functional has been implemented in 
the linear-response approximation (8,9) of Vignale and 
Rasolt [7, 8]. Three different parametrizations for the 
orbital susceptibility entering this functional have been 
tested. The LCH 3-term parametrization was found to 
be slightly less accurate than claimed in the original ref- 
erence (rms error 1.9 x 10~ 3 instead of 1.5 x 10~ 3 ). Two 
alternative parametrizations were developed. Our 3-term 
parametrization leads to an rms error of 1.2 x 1CP 3 , while 
the 5-term expression we used in our numerical calcula- 
tions reduces this to 2.1 x 1CP 4 . 

(iii) Orbital currents in open shells result in small but 
not negligible shifts of the ionization energies. Such shifts 
can become important, e.g., when ionization takes place 
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in the presence of external electric or magnetic fields. The 
same applies to the formation and breaking of chemical 
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calculation of excitation energies. The calculated shifts, 
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over experimentally. A more conclusive test of (p)CDFT 
would thus require experiments carried out in the pres- 
ence of magnetic fields, giving rise to a well defined cur- 
rent in the open shell. 

(iv) The order of magnitude of the pCDFT terms 
shows that modern density-functional techniques have 
reached an accuracy at which the magnitude of small 
current-related effects arising in open shells is beginning 
to be significant compared to the remaining difference 
to experiment. Hence further refinement of DFT-based 
calculations of atomic spectra should consider the possi- 
bility of spontaneous or induced currents in open shells 
[17, 19, 25, 30]. 

More research is needed to extend this analysis to other 
atoms and to molecules. Fully self-consistent (neither 
perturbative nor post-LDA) CDFT calculations would 
be desirable to this end, as is, obviously, construction of 
more reliable current-dependent xc functionals. 
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